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Abstract 

We investigate the backreaction equations for an acoustic black hole formed in a 
Laval nozzle under the assumption that the motion of the fluid is one-dimensional. 
The solution in the near-horizon region shows that as phonons are (thermally) radi- 
ated the sonic horizon shrinks and the temperature decreases. This contrasts with the 
behaviour of Schwarzschild black holes, and is similar to what happens in the evap- 
oration of (near-extremal) Reissner- Nordstrom black holes (i.e., infinite evaporation 
time). Finally, by appropriate boundary conditions the solution is extended in both 
the asymptotic regions of the nozzle. 
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1 Introduction 



In 1974 Hawking [1] announced to the physical community his famous result: black 
holes are not black at all, but emit thermal radiation at a temperature proportional 
to the horizon surface gravity. Hawking came to this remarkable conclusion by exam- 
ining, within the framework of quantum field theory in curved space, the propagation 
of a massless scalar field in the geometry of a collapsing body forming a black hole. If 
one assumes spherical symmetry, the space-time in the region exterior to the body is 
described by the Schwarzschild solution (Birkhoff theorem). For these kinds of black 
holes the horizon is located at r = 2M (M is the black hole mass) and its surface 
gravity is k — 1/4M. This yields an emission temperature T H — || = (the 
velocity of light and Boltzman constant are set to one). As a consequence black holes, 
the final end-state of gravitational collapse according to classical General Relativity, 
are unstable. 

The actual evolution of an evaporating black hole should in priciple be described (un- 
til quantum gravity effects become relevant) by the semiclassical Einstein equations 
(backreaction equations) 

G/xAgap) = 8vr (T^^)) , (1.1) 

where G^ u is the Einstein tensor and the r.h.s. is the expectation value of the quantum 
stress tensor for the matter fields which drive the evaporation. 

Eq. (1.1) has to be solved for the evaporating black hole metric g a p. This requires 
the knowledge of {T^ y {g a p)) for an arbitrary (let's say spherically symmetric) metric. 
Unfortunately, because of the extreme difficulty of the problem, no such expression 
is available and the backreaction is usually modelled by extrapolating Hawkings's 
result, which is strictly valid only for static or stationary black holes. Therefore a 
spherically symmetric black hole of mass M(t) is supposed to emit radiation at a 
temperature T = h(8nM (t)) -1 . Because of this emission, it looses mass at a rate 
given by the Stefan-Boltzman law ^ ~ —AT 4 ~ -1/M 2 , where A = IQttM 2 is the 
horizon area. As the mass decreases the black hole becomes hotter, as can be seen 
from the temperature formula given above, and eventually disappears (with a final 
explosion?). Its lifetime is of the order Mq where Mo is the initial mass of the black 
hole. 

If the black hole possesses a conserved electric charge Q (\Q\ < M) the picture is 
quite different. In this case the space-time is described by the Reissner-Nordstrom 
solution for which the surface gravity reads 

K = (r+- M)jr\ (1.2) 

where r + is the radius of the horizon, i.e. r + = M + \fM' 2 — Q 2 . One sees that, for 
near-extremal black holes (M > \Q\), as the mass of the hole decreases because of 
the evaporation, the temperature T = -^k drops costantly. When M is reduced to 
M — \Q\ the temperature is zero. This final state (remnant) is reached in infinite 
time. This corresponds to the so called third law of black hole thermodynamics. 
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In 1981 a remarkable paper of Unruh [2] appeared showing that a quantum emis- 
sion similar to the one predicted by Hawking for black holes is expected in a seemingly 
completely different physical context, namely fluids undergoing hypersonic motion. 
This far reaching result opened a continuously developing field of research (the so 
called black hole analogue models [3]) in condensed matter physics where, unike grav- 
ity, the hope to perform experimental test on these theoretical predictions does not 
seem so remote. 

In this paper we will give a first insight on the backreaction this emitted radiation 
has on the fluid. The plan of the paper is the following. 

In section II we report Unruh's analysis of the sound propagating in a hypersonic fluid 
within the action formalism which is then used to write the backreaction equations 
for the fluid motion. 

In section III we introduce a lower dimensional model where the backreaction equa- 
tions can be solved to find the first order in % correction to the classical flow. 
In section IV we evaluate the quantum stress tensor for the classical background to 
be inserted in the backreaction equations. 

In section V we solve the backreaction equations near the sonic horizon. 

In section VI and VII the near horizon solution is extended in the asymptotic regions 

through the imposition of appropriate matching conditions. 

Section VIII is devoted to the conclusions and some extrapolations of our results. 
In Appendix A and B we report tecnical details needed to perform the calculations. 



2 Hypersonic flow and sonic black holes 



Unruh considered an irrotational, homentropic fluid. In this case the Eulerian equa- 
tions of motion can be derived from the action 



S 



(fx 



+ -p(W>) 2 + u{ P ) 



(2.1) 



where p is the mass density, the velocity potential, i.e. if = and u the internal 
energy density. The overdot means differentiation with respect to (newtonian) time. 



Varying S with respect to ip yields the continuity equation 

p + V ■ {pit) = , 
whereas variation with respect to p gives Bernoulli's equation 

ij + ^ir 2 + p{p) = o , 



(2.2) 



(2.3) 



where p(p) = %. 
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One can now obtain the linearized wave equation for the propagation of sound 
waves in a background mean flow by replacing 



P 



P + Pi , 



(2.4) 



where p and ip define the mean flow; they are assumed to satisfy the equations of 
motion and ipi and p\ are small amplitudes perturbations. 
Expanding S up to quadratic order in these perturbations one obtains 



S = S - j d A x 

= So + S 2 . 

The speed of sound c is defined as 



1 c 1 

piti + ^-p\ + ^WO 2 + prf ■ V^i 



(2.5) 



2 djJL 

^ dp 



(2.6) 



One can derive the equation of motion for p x 

= -4(^i+^-^i) • 



pi 



(2.7) 



Since p\ occurs quadratically in (2.5) we may use eq. (2.7) to eliminate it and obtain 
an action for the potential ifj\ only 



S 2 = -J d*x [ip(^) 2 - + V • vVi) 2 

The remarkable thing is that S 2 can be written as 

S2 = ~\S ^v^fT^i^i , 

where, following Unruh, we have introduced the acoustic metric 

„2 



9nu 



pi c z - v 2 lf T 
c\ if -I 



(2.8) 



(2.9) 



C2-10) 



and / is the three-dimensional identity matrix. 

In the form of eq. (2.9) S 2 is completely equivalent to the action for a massless scalar 
field ipi propagating in a curved space-time whose line element is 



ds 2 = g^dx^dx" 



(2.11) 



Varying the action (2.9) with respect to ipi gives the sound wave equation which can 
be written in a simple and elegant geometrical way 



□^i = 0, 
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Figure I: A Laval nozzle. The waist of the nozzle represents the sonic horizon (\v\ = c). 
In the region on the right of the waist \v\ < c and on the left \v\ > c (sonic hole). 



where □ = V M V M is the covariant derivative with respect to the acoustic metric g^ u . 
Variation of S 2 of eq. (2.9) with respect to the acoustic metric gives the so called 
"pseudo energy momentum tensor" (PEMT) T M „ 

r - = -^- (2 ' 12) 

As discussed by Stone [4] the invariance of S 2 under diffeomorphisms implies the 
covariant conservation of T^ v 

T^ ;iy = 0, (2.13) 

which describes in a nice compact form the exchange of energy and momentum be- 
tween the waves and the fluid reproducing various known laws of classical fluid motion. 
Of particular interest are fluids undergoing hypersonic motion. In this case the region 
of the fluid for which if 2 > c 2 is called acoustic black hole, its boundary if 2 = c 2 
defines the sonic horizon. From this region sound waves cannot propagate upstream 
since the fluid velocity is bigger than the sound velocity. Sound is dragged by the 
fluid and cannot excape. 

Using the same arguments of Hawking, Unruh, quantizing the field ipi, showed that 
in the formation of a sonic hole one expects a thermal emission of phonons at a 
temperature 

T= — , (2.14) 
where k is the surface gravity of the sonic horizon [5] 



(2.15) 

H 



n is the normal to the horizon. 



A tipical situation where this is supposed to occur is a Laval nozzle depicted in 
Fig.l. The fluid flows from right to left; at the waist of the nozzle the fluid velocity 
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reaches the speed of sound. This is the location of the sonic horizon. 

Expanding into the left region the velocity increases and the fluid becomes hypersonic; 

this region describes the acoustic black hole. In the asymptotic right region a thermal 

radiation of phonons is expected according to Unruh's analysis. 

The aim of this paper is to give a first qualitative description of the effects this emitted 

radiation has on the fluid dynamics, i.e. the backreaction for acoustic black holes. 

To this end we have to establish the analogue of the semiclassical Einstein equations 

(1.1), i.e. the evolution equation for the fluid driven by the quantum fields. 

Using standard background field formalism one writes the fundamental quantum fields 

as a sum of background fields (p, ip) (not necessarily satisfying the classical equations 

of motion) plus quantum fluctuations. Integrating out the quantum fluctuations one 

obtains the one-loop effective action formally defined as 

r = S + hiTr In [n gM ] +0(h 2 ) , (2.16) 

where S is the classical action (2.1) and n g ^ p ^ is the DAlembertian calculated from 
an acoustic metric g^p^v) of the functional form of eq. (2.10) [6]. We assume 
the divergences in the determinant of the above effective action to be removed by 
employing a covariant regularization scheme. This is the key hypothesis of our work. 
We will comment on it in the conclusion section. One can therefore write r = S + 
S q (g^u(p,v)), where the quantum part of T depends on the dynamical variables p, v 
only through the acoustic metric (this is a direct consequence of our hypothesis) and 
coincides with the effective action of a massless scalar field propagating in a curved 
background whose metric is g^v(p,v). 
Given this it follows that, using the chain rule, 

6S, _ 5S ± 5jT__VEl {T ) S JH-VZI (T (2 17) 

5p ~ 6<r Sp - 2 Sp - 2 {1 ^{p) 9 [Z - U) 

2 p 2 c 



where 



and (T) = (Tp. 
Similarly 



SS q _ 6S q 6<r ^~g 6<r (9 1Q , 



5vi 5g» v 5vi 2 Sv 

= ^(T tl )+v t (T lt )) . 

In the above expression g^ v is the inverse metric, i.e 




(2.20) 
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The backreaction equations are therefore (taking into account that does not de- 
pend on ip) 



p + ^ • {pt) - V, (f {(T ti ) + v t (T u ))) = 
i> + \-t 2 + P(p) - H(T> = . 



(2.21) 
(2.22) 



The (T^v) appearing in the backreaction equations should be considered as expectation 
value of the quantum pseudo energy momentum tensor operator taken in a quantum 
state appropriate for the radiation process, namely the analogue of the Unruh state 
[7] (the quantum state suitable to describe black hole evaporation at late times), in 
which in the remote past prior to the formation of the sonic hole the quantum field is 
in its vacuum state. 

Inspection of the backreaction equations reveals that the continuity equation (asso- 
ciated to the symmetry — > + const, of the effective action) gets modified by 
including the divergence of the phonons momentum density (y/—g(T®)), whereas the 
Bernoulli equation gets an additional contribution to the chemical potential given by a 
term proportional to the trace of the pseudo energy momentum tensor of the phonons. 
These equations reflect the underlying two-component structure of the system as in 
the Landau-Khalatnikov theory of superfluidity [8]. 

As in the black hole case the general expression of S q and (T^ v (g a p)) is completely 
unknown and no explicit solution of the backreaction equations (2.21), (2.22) can be 
given. 

3 The dimensional reduction model 

Lower dimensional models are very useful in physics as they allow explicit solutions 
of the dynamical equations to be found, providing at least a qualitative description of 
the evolution of real four-dimensional systems. This attitude has been largely used in 
studying quantum effects in black hole physics, where two-dimensional models enable 
investigation of the so called "s-wave sector" [9]. The most popular one is the one 
proposed by Callan-Giddings-Harvey-Strominger (CGHS) [10] in which the 4D quan- 
tum stress tensor is replaced by a 2D one associated to a minimally coupled massless 
scalar field described at the quantum level by the Polyakov action. This approxi- 
mation neglects backscattering of the quantum fields caused by the potential barrier 
located near the horizon. 

Within this philosophy a qualitative insight in the backreaction going on in a hy- 
personic fluid can be obtained by assuming a one dimensional flow for the fluid; the 
relevant physical quantities depend only on the z coordinate running along the axis 
of the Laval nozzle of Fig. 1 (i.e. velocity components along x and y are negligible 
with respect to the z component). 
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The effective action for the CGHS-like model for the fluid quantum dynamics can be 
given as r^ 2 ^ = S^ 2 ) + S po i, where 



d 2 xA 



pi) + -p{d z %l)f + u{p) 



(3.1) 



is obtained integrating S over the transverse coordinates x, y and A is the area of 
the transverse section of the nozzle. Using the chain rule as before one can write the 
backreaction equations as 



Ap + d z (Apv) + d z 



--((T{?)+v(T^)) 







^ + y + Mp))4(< T(2) >)= ' 



where 



<7i 2) > 



SS, 



pol 



(2) 

and g ab is the (t, z) section of the acoustic metric (2.9) 



a {2) 

9ab 



P ( C 2 —v 2 V 
C \ V -1 



(3.2) 
(3.3) 

(3.4) 



(3.5) 



Our aim is to solve these backreaction equations to first order in h to calculate 
the first quantum corrections to the classical fluid configuration. Therefore we need 
the quantum source terms calculated at zero order in H, i.e. for the classical solution. 



4 The quantum pseudo stress tensor for static sonic 
black holes 

Before undertaking the study of the backreaction in eqs. (3.2), (3.3) let us describe 
the solution of the classical equations of motion which describes a static acoustic black 
hole and then find the corresponding (T a5 ) in the Unruh state for this background. 
The classical equations of motion are given by eqs. (3.2), (3.3) neglecting the phonon 
quantum contribution. 

Assuming steady flow, i.e. no time dependence, the continuity equation (3.2) yields 

A(z)p(z)v(z) = const = D , (4.1) 
whereas the Bernoulli equation simplifies to 

y + Mp)=0. (4.2) 
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Differentiating both equations and taking into account the definition of the velocity 
of sound (eq. (2.6)) one easily gets 



A „. V X J?. ) ' 



(4.3) 



which is the famous nozzle equation. From this one deduces important features of 
the fluid motion. First of all the horizon (\v\ = c) forms at the waist of the nozzle 
(A' = 0).* Furthermore in the hypersonic region (\v\ > c) the velocity increases to- 
gether with the section of the nozzle, just the opposite of what happens in the more 
familiare subsonic (|i>| < c) case. 

Let us further assume the sound velocity c constant. This allows the Bernoulli equa- 
tion (4.2) to be written as 

2 po 



0, 



i.e. 



p = p e ^ . 
Let us choose the profile of the Laval nozzle as 

A= [ A + f3z 2 ior\z\ < z 
1 A + = const for | z \ > z 

where A , (3 and z are constants. 

From the continuity equation (4.1) one has for \z\ < z 



(4.4) 
(4.5) 



(4.6) 



z — — 



D e v2 ' 2c2 
Po v 



-Ac 



(4.7) 



The sonic horizon is located at z = where v — — c (remember that the fluid is 
moving from right to left). The constant D in eq. (4.1) can be evaluated as 



-A p e 2 C = D , 
which when inserted in eq. (4.7) yields 



(4.8) 



(4.9) 



This implicitly defines v = v(z) and from eq. (4.5) we have p(z); so for \z\ < z$ the 
acoustic metric g^yz) of the static sonic black hole is determined. 
For \z\ > z the solution is simply v(z) = const = v(z ) for z > z and v(z) = const = 
v(— zq) for z < —zq. Similarly for p(z). Note that \v(z)\ increases as z decreases, while 

*Note that this is in general not true if external forces are present. In that case, the equation 
giving the location of the sonic horizon is A' /A = —f/pc 2 , where / is the external force (see the 
second of Refs. [14]). 
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p{z) decreases as z decreases. This resulting acoustic metric should be regarded as the 
asymptotic configuration of the fluid resulting from the (time dependent) formation 
of the sonic hole (see for example [11]). 

This dynamical evolution of the fluid excites the modes of the quantum field, which 
is supposed to be in its vacuum state before the process of formation of the sonic hole 
begins, yielding nontrivial expectation values for the quantum PEMT. At late times 
these may be approximated by (T^) in the Unruh state evaluated in this background. 



We introduce a set of null coordinates 

x~ = c 1 1 - 



X 



cit + 



dz 
c + v / 
dz 



c — v 

In terms of these coordinates the acoustic metric is 

p c 2 - v 2 ( \ 



a {2) 

Sab 



c cr 



2 



with a, b = x + , x_ and the quantum part in eq. (3.2) reads 



— ((T tz )+v(T zz )) = -c 2 
c 

The Polyakov (T ) in the Unruh state is [12] 

% 



(T-) + (T ++ ) 



(c + v) 2 (c — v) 2 



(T__> 



12tt 



48tt 



whereas 

(T ++ ) - 

Here the conformal factor of the metric C is 



12tt ' ++ 



= P c 2 - v 2 
c c 2 



and k is related to the surface gravity k on the horizon 



d 



k = C K = c—v 
dz 



H 



For our particular chioce of A(z) (see eq.(4.6)), the surface gravity is k — 
c 2 \J f3/A Q . The trace (T^) is completely anomalous and is given by 

/ T (2)\ JLd(2) 

K ; 24tt 



(4.10) 
(4.11) 

(4.12) 

(4.13) 

(4.14) 
(4.15) 

(4.16) 
(4.17) 



C 2 K = 



(4.18) 



9 



where is the Ricci scalar for the metric gf b . 

Note also that in the Unruh state (T ) is regular on the horizon (where v = —c) 

making the r.h.s. of eq. (4.13) finite. 

Explicit calculation of the components of (T^) are given in Appendix I. We recall 
here only the asymptotic limits 

(T-)U 
{T ++ )\ z - 



48tt ' 

flK 2 

~48^ 



(4.19) 



5 The backreaction equations near the horizon 



As said before our aim is to find the first order correction in % to the classical sonic 
black hole solution p(z),v(z) described in eqs. (4.5)-(4.7). To this end we define the 
quantum corrected velocity potential and density by 



Pb = p(z) +ep 1 (z,t) , 



(5.1) 
(5.2) 



where the subscript "B" stays for backreaction, v = d z ip and we have introduced a 
dimensionless expansion parameter e = h / (\D\A ) . These values are found by solving 
the backreaction eqs. (2.21), (2.22) linearized in e 



e [Ah + d z (A( Pl v + pdM)] = c 2 d z 



+ 



(T. 



++/ 



(c + v) 2 (c — V ) 



+ vd z ijj 1 + —pi) 



= eGo(z) 



tF 2 {z) , (5.3) 
(5.4) 



The explicit expressions for F 2 and G 2 are given in Appendix I. 
From eq. (5.4) one can obtain p\ 



P 



— ~ "01 - «Vl 



(5.5) 



where we indicate ' = d z . 

Taking the time derivative of this equation and substituting it in the continuity eq. 
(5.3) one obtains the following wave equation for ipi 



= H{z) . 



(5.6) 



The l.h.s. is simply proportional to V^^^V^V'i) and the r.h.s. H(z) represents 
the quantum source evaluated for the classical static solution. 
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We are mainly interested on the fate of the sonic horizon, so we will solve the wave 
equation (5.6) near z — 0. To this end we need an expansion of the background 

(2) 

quantities p and v, and hence {T^ b ), for kz <^ 1. The results of this tedious calculation 
are given in Appendix II. We give here only the final results. Eq. (5.6) becomes 

+ KZ)^! - 2lp[ - 2c«(l + §K2)^i - = 

= H + f^x/cz , (5.7) 

where 

3 2/1 

# ~ 0.038^^ , (5.8) 

7T 

Hl ~ 1.128^^. (5.9) 

7T 

We shall write the solution in the following form 

2 2 

V>i = / (t) + Kzh(t) + ^-/ 2 (f) + ... (5.10) 



where fi(t),i = 1,2,3... are functions of t only, to be determined. 
The zero order (in kz) equation is 

-J (t) - 2«A(*) - 2c/€ 2 /i(t) = #o , (5.11) 



whereas the first order equation is 

-\fo(t)+f 1 (t)-2Kf 2 (t) 



Ack 2 



h(t) + lh(t) 



Hi . (5.12) 



We shall now define the boundary conditions on the solution. We require that at 
some given time, let's say t — 0, the evaporation is switched on, i.e. we require 
fa(t = 0) = ui(t = 0) = pi(t = 0) = 0. 

The vanishing of ip\ and vi = d z ipi at t — implies from eq. (5.10) that = 0) = 0. 
Remembering now the relation between p\ and ip\ (see eq. (5.5)) we have 

Pl (t = 0) = ^(t = 0) = ^ i.e. 

where the expression on the r.h.s. means the coefficients of the i-th term of the 
expansion in z of G 2 j A. 

Evaluating eqs. (5.11), (5.12) and the time derivative of eq. (5.11) at t — one can 
obtain the following approximate expression for ipi (see Appendix II for details) 

7 2 2 

QjO n Of! o Oon .HZ, . .. 

^1 = Oit + yt 2 + ^yt 3 + .... + Kz{ht + yt 2 + ....) + —(Cit + ...) . (5.14) 
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The numerical coefficients are reported in the same Appendix. 

Since we have calculated the quantum source only for the static classical metric 
there is no time evolution of the source itself. This limits the validity of our solution for 
small values of t only (cut < 1). For this reason we have just given a power expansion 
of the functions fi(t). Therefore we are only able to predict how the backreaction 
starts, i.e. to connect the static fluid configuration at time to = to a quasi static 
fluid configuration at time to + At. With this in mind we can now give the quantum 
correction to the velocity field 

Vi(z, t) = Kbit + (nz)ciKt = Kt(b\ + C\kz) . (5.15) 

Note that b\ > (see Appendix II) and being kz <C 1 for the validity of our solution, 
we can conclude that V\ > 0. The background velocity is negative, so we conclude 
that the backreaction has the net effect to decrease the modulus of the velocity, i.e. 
the fluid is slowing down. 
The quantum corrected velocity is therefore 

vb = v + evi = — c + ckz ck 2 z 2 + e(bi + C\Kz)Kt . (5.16) 



The acoustic horizon is defined \vb\ = c; this yields 

CK + hCiK 2 t C 

This equation shows that the horizon is moving to the left with respect the classical 
location zh = 0. Therefore, as the evaporation proceeds, the hypersonic region gets 
smaller and smaller. This is the behaviour one would have naively expected. The 
coefficient b± determining the quantum correction to the velocity and hence the evolu- 
tion of the horizon is just the gradient of the additional chemical potential related to 
the expectation value of the trace evaluated at z — 0. This should be compared to the 
black hole case where the evolution of the horizon is determined by the energy flux 
(M oc (T r t ) in spherical symmetry). While in the latter case Hawking radiation oc- 
curs at the expense of the gravitational energy of the black hole, in the fluid phonons 
emission takes away kinetic energy from the system. 
Finally, from eq. (5.5) we can calculate p\ 

_ i 

pi = ^—z — (a + 5nz)t , (5.18) 

where 

a = —a 2 + CKb\ , (5.19) 
S = —b 2 — a 2 + ckc\ . (5.20) 



This yields p 1 < 0, i.e. Pb = P + e Pi is decreasing. 
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6 The backreaction equation in the region z > Zq 



In the region z > zq, assuming kzq <C 1, the background solutions can be approxi- 
mated as 

V — — C + CKZq , (6.1) 

p = p e~^(l + kz ) . (6.2) 

The backreaction equation (5.6) becomes the simple flat space wave equation 

d_d + ipi = (6.3) 

Since in this region (T^) = (T+_) = and (T__) = const, then F 2 = G 2 = 0. 
The null coordinates x~,x + are defined as follows (see eqs. (4.10), (4.11)) 

x~ 

x + 

The most general solution of (6.3) is a linear combination of two arbitrary functions 
of the kind £(x~) and i](x + ). 

The source term given by the quantum pseudo stress tensor is discontinuous in 
z = z , so our model resembles a sandwich of space-time regions glued across singular 
hypersurfaces. We can find a solution of (6.3) for z > z by requiring continuity for 
the vi(z,t) field across the boundary z = z . Because of the discontinuity in (T^ b ), 
there is no way to make v\ and p\ simultaneously continuous across z . Our choice 
leads to a quantum corrected acoustic metric wich is continuous in z = zq up to a 
conformal factor. 
For z < zq we have 

Vi(z, t) = {pi + C\KZ)Kt (6.6) 

and therefore 

Vi(zQ,t) = (bi +CiKZo)Kt . (6.7) 

As seen from the curved space point of view, the z = z surface is a timelike surface 
of the acoustic metric of eq. (3.5). It is crossed in the future only by outgoing rays. 
Therefore for z > zq the solution of the wave equation assumes a retarded "Vaidya 
form" ipi = ipi(x~). Now from eq. (6.4) 

vi = d^xipT) = — l —d^ x = F(x~) . (6.8) 

KZq 

Taking into account that x~(t = 0, z = zq) = —1/k = Xq, F(x~) can be approximated 

as 

F(x~) = l(x~ - Xq)Q(x- - xq ) + 0((x~ - xq) 2 ) , (6.9) 



c[t- 



z 



■ , 



c[t + 



CKZqJ 

z 



2c — CKZq 



(6.4) 
(6.5) 



13 



where / is a constant to be fixed by requiring continuity in z = z : 

F(x~) = Id = Vi{z Q ,t) = {pi + CKZq ) Kt , (6.10) 

Z = ZQ 

so that 

I = —(pi + CiKZq) , (6.11) 

c 

which is positive definite. 

Therefore for z > z , x~ > Xq , we have 

vi = vi(x~) = ~( K h + c 1 k 2 z )(x' - Xq) + ... , (6-12) 

whereas for x~ < Xq it vanishes identically. 

The quantum corrected velocity field reads therefore, for z > zq, 

vb = v + hvi = — c(l — kzq) H — (bi + C\KZq)(x~ — Xq) , (6.13) 

c 

showing a velocity decreasing in modulus with the advanced time x~ . 
Integrating vi = d z ipi(x~) one can get ipi(x~) 

^1 = --7r( Kh l + C 1 K 2 Z )(X~ - Xq f (6.14) 

and from it 

_ i 

Pi = ~ VV t) = P °^2 ^ + KZ 0)( K bl + CiK 2 Z )(x" - Xq )Q(x~ - Xq) , (6.15) 

which increases in terms of the advanced time. 

As already said p\ cannot be made continuous across z = z . Evaluating the limits 
from both sides of the hypersurface z — z one can find the jump in p 1 

_ i 

Api = — [a 2 + KZQ^bx + b 2 + a 2 )ct] . (6.16) 

cr 



7 The backreaction equation in z < —Zo 

In this region the backreaction equation is d-d+tpi = 0, as in the previous case. Now 
the null coordinates are 

x~ = c(t + —) , (7.1) 



CKZq 

2c + CKZq 



* + = C [ t + w~~~zr ) • (7-2) 
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As one can see from the acoustic metric (3.5), the surface z = —z is now a spacelike 
surface; it is crossed in the future direction from both ingoing x + = const and outgoing 
x~ = const characteristics of the wave equation. So the solution will depend on both 
x~ and x + . 



The classical background quantities are now for z < —z 

V = — C — CKZq , 
_ 1 

p = p e 2(1 - K z ) . 



(7.3) 
(7.4) 



The boundary conditions we will impose to find the solution in this region are the 
continuity for both v\ and p\. Unlike the z = z case, here this is possible being 
z = —z a spacelike surface. Now we impose continuity of t>i 



d z ipi(-z ) = <9 2 Vi(-^) = <Pi(ct) , 



(7.5) 



where the r.h.s. is calculated on the solution obtained in Section 5. Continuity of pi 
reads 



P 



— [-Ip! - Vd z lj}i] 



(7.6) 



Being the background quantities (7.3), (7.4) continuous across z = —z and taking 
account of eq. (7.5) we require 



d z ^i(-z ) = h/>i 



cip 2 (ct) . 



From Section 5 



(fl(ct) = -(&! - C X KZ^)t , 
C 

/ ,x fa 2 b 2 KZ \ 
ap 2 {ct) = ct 



(7.7) 

(7.8) 
(7.9) 



and using standard technique for solving the wave equation one eventually arrives at 

(x~ -x )Q{x- -Xq) 



Vl 



-^-(bl ~ CxKZq) - ^(^2 - b 2 KZ ) 



-(x + - x^)Q(x + - x^) + -(&i - c 1 kz )(x -x )Q(x - z ) , (7.10) 



where we have defined 



X r 



x (z — —z , t = 0) = — 



K 



Xq = X + (z = — Zq, t = 0) = 



2 + kzq 



(7.11) 
(7.12) 



15 



So we have a Vaidya ingoing solution for x < x 



Vi = - 
whence for x" > Xq 



KZq ( j 2 1 



. («6i - 2 K Ci) - ^2"( a 2 - «2o&2) 



- x+)6(a; + - x+) , (7.13) 



«i = 



KZq 



(nbi - CiK z Q ) -(a 2 - b 2 KZ ) 



H (61 - C 1 KZ )(x' - Xq) 

c 



z + z 



KZ (2 + KZ ) 



In a similar way we obtain 



Pi 



(bi — ciKZ )(x — x )Q(x - x ) + 



K 2 Z 



(6l - CxKZq)- 



-(a 2 - b 2 KZ ) 



[x — x )&(x — x ) — (x — Xq)Q(x — X, 



»]}■ 



(7.14) 



(7.15) 



8 The fate of the acoustic black hole 



The aim of this paper was to face the backreaction problem for acoustic black holes. 
First we have found the equations for the fluid driven by the linear quantum fluctua- 
tions, i.e. the linearized "backreaction equations" . 

In a lower dimensional case, keeping the Polyakov approximation for the stress tensor, 
we have solved these equations near the horizon, under the assumption that the quan- 
tum source corresponds to the pseudo stress tensor of the quantum fields evaluated 
for the static classical background. 

The information we can gain is just preliminary: we obtain an indication on how the 
quantum solution starts to depart from the classical configuration, once the quantum 
effects are switched on. 

It's not possible to do anything better even for the gravitational black hole evapora- 
tion. However very interesting indications seem to emerge even from this simplified 
model. 

We have shown that the sonic horizon moves to the left with respect to the nozzle and 
the hypersonic region shrinks. From the expression of the quantum corrected fluid 
velocity 

v b = — c + ckz — -ck 2 z 2 + e(bi + c\Kz)nt (8.1) 



6 



and the equation for the horizon 



zh 



— e — t 

c 



(8.2) 
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one can evaluate how the emission temperature of these quasi static configuration 
varies with time, i.e. 



T = 



he 8vb 



2tt dz 

he ( 
" 2^ I 



1 2 €K 2 Ci N 

K - -K Z H H 1 

6 C , 



(8.3) 



This, using eq.(8.2), yields 

T 



he 
— / 

2vr 



-K 



1 + 



he 
2tt 



tK ( bi 
~~c \~3 
563 



720tt 



+ cij t 

€K 3 cA t 



(8.4) 



This expression is rather significative: it tells us that the emission temperature of the 
sonic black hole decreases in time. Unlike Schwarzschild black holes, the sonic black 
hole gets cooler as it evaporates. Its behaviour is therefore much more similar to a 
Reissner-Nordstrom black hole. 

It's very interesting to try to extrapolate this behaviour in time. From eq. (8.4) 
we have 

AT oc -n 4 c 2 A At ~oc -T 4 At . (8.5) 
This yields approximatively 



This implies that as the sonic black hole radiates, its temperature decreases reaching 
asymptotically zero value in an infinite time. 

For a near extremal Reissner-Nordstrom black hole (r + ~ M ~ \Q\) a similar 
analysis yields 



dM 



dt 



M 



~ —AhTJj oc —Q (M — Q) Q , 
1 

— = t ~ — . 
Q T 2 



(8.7) 
(8.8) 



So sonic balck holes resemble near extremal R-N black holes. The end-state of the 
evaporation process corresponds in both cases to a zero temperature configuration, 
reached asymptotically in time. 

One should mention that other analog models like a thin film of 3 He- A with a moving 
domain wall [13] seem to show a non-vanishing end-temperature of the evaporation 
process. 

In a similar way one can treat sonic black holes formed by a Bose-Einstein condensate 
[14]. The major difference comes from the fact that the sound velocity c is not 
constant, but proportional to pa. This will be discussed elsewhere. 
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9 Conclusion 



The basic question which the analysis presented here leaves unanswered is to what 
extent does the behaviour we have obtained for the quantum corrected evolution 
depend on the linear dispersion relation used (free scalar field to describe phonons), 
which ignores short distance corrections due to the molecular structure of the fluid. 
The same uncertainty exists for any study in curved space quantum field theory which 
involves high frequency modes like the Hawking effect or in general renormalization 
calculations. One does not know how the yet unknown behaviour of the spacetime 
at the Planck scale may influence the results. So far many studies have been devoted 
to investigate the robustness of the thermal radiation predicted by Hawking both 
in the gravitational context and in the analogue models. Looking at 2D models 
one has mainly analyzed the behaviour of the Bogolubov f3 w coefficient giving the 
created particles number once nonlinear dispersion relations are introduced in the 
high frequency regime. This is supposed to model the breaking of local Lorentz (or 
Galilean, in the analogue model context) invariance here expected. The scale at 
which this happens is the Planck scale for gravity, whereas for sonic black holes it 
is the intermolecular size. As a result of these studies there is now some consensus 
(see however the recent paper [15] where interesting counterexamples are given) that 
the (3 W coefficient is basically unaffected by the dispersion relation [16]. However no 
one knows yet what happens to observables like (T^), which are in addition sensible 
to the regularization scheme used. In a pure 2D context and for massless scalar 
fields, Jacobson [17] has argued that, within a covariant regularization, no significant 
deviation from the usual expression for the trace anomaly (eq. (4.18)) and the flux 
(eq. 4.19)) are expected if one introduces a cutoff at high frequencies. However, 
for the hydrodynamical system we have considered, covariance is a symmetry of the 
phonons low energy effective theory only, which is broken at short distance. Hence non 
covariant terms depending on the microscopic physics are expected to show up in the 
effective action and are crucial for a correct description of the unperturbed quantum 
vacuum of the fluid. However the expectation values (T^) entering the backreaction 
equations (2.21, 2.22) do not represent the energy momentum of the fluid quantum 
vacuum. They describe instead the perturbation of the stationary vacuum (whose 
energy is strictly zero [18]) induced by inhomogeneities and by the time dependent 
formation of the sonic hole which triggers the phonons emission. In this paper we 
have assumed that these deviations can be computed within the low energy theory. 
This situation is not unusual. Casimir effects are well known examples of vacuum 
disturbances caused by the presence of boundaries. It happens that the Casimir energy 
is often (but not always, see G. Volovik in Ref. [3]) independent on the microscopic 
physics and can be calculated within the framework of the low energy theory. This 
happens because, while low frequency modes are reflected by the boundaries, for the 
high energy ones the wall is transparent. They produce a divergent contribution to 
the vacuum energy which is canceled by a proper regularization scheme and does not 
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affect the finite result. We have assumed that a similar decoupling happens for the 
acoustic black hole. The check of our hypothesis would require an analysis of the 
quantum system within the microscopic theory which takes into account the time 
dependent non homogeneous formation of the sonic hole. This is for the moment 
beyond computational capability. 

More efforts are therefore necessary to fully understand the role played by the high- 
frequency modes in the above scenarios; in particular what is needed is a proper 
description of the way these modes interact with the underlying medium. 
To this end the study of black hole analogue models seems much more promising, 
since there the underlying physics is at least in principle known. Progress in this area 
may then be used as "theoretical laboratory" for ideas to be exported and tested in 
quantum gravity. 

Given the the present limited understanding, our work is of course not the end story 
but it should be regarded simply as a first attempt to attack the backreaction problem 
in sonic black holes, an attempt performed using available technical tools developed 
in quantum field theory in curved space. 

It is obvious that more detailed analysis are needed and will hopefully be performed 
in the not too distant future. However being the behaviour we have found for the 
evolution of a sonic hole physically quite reasonable, we would expect these features 
to be confirmed. 
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In this appendix we give the relevant expressions needed to evaluate the r.h.s. of the 
backreaction equations. 

The stress tensor components are 



A Appendix I 




n (c 2 -v 2 ) 

127T 4C 4 



( 




2p 4 p 2 



(c 2 - v 2 ) + 




(1.1) 




Hk 2 



(1.2) 



++ 



48tt 
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The anomalous trace is given by 



where 



R = — 



pc 



T (2) = AflP) , 

24vr 



;i.3) 



(1.4) 



One needs an expansion of the background velocity field till the fifth order in kz 



1 9 9 11 o o TT 

V ~ — C + CKZ CK Z CKT H 

6 36 1080 



4 4 , 



769 
4320' 



CK 5 Z 5 



Similarly for the energy density 

_i / 1 2 2 23 3 3 167 4 4 485 5 5 

p - p ° e 2 [ 1 + KZ - - 6 ^ 2 - ^ + Ym + Wa k z 



(1.5) 



;i.6) 



Inserting the expansions (1.5) and (1.6) it is possible to obtain the following expres- 
sions, entering in the backreaction equations 



G 2 ~ 
H{z) ~ 



K 2 C 2 A 2 / 




24tt V 


5 + 


|£>|A)« 3 


( J 


48vr 




k 3 c 2 A / 


109 


48tt V 


60 



49 



2 5 
649 4909 



2 2 
K Z 



60 360 

19501 
+ ^6T^ 



-KZ 



(1.7) 
(1.8) 
(1.9) 



B Appendix II 



The coefficients entering in the solution of the backreaction equations are given as 
follows 



Oi = 


/o 


a 2 = 


/o 


a 3 = 


/o 


6i = 


A 


&2 = 


A 


Cl = 


A 



t=o 3 7 
1189 



t=o 



t=o 



t=0 



-'yen 
t=o 120 

10337 



360 



7 (c/t) 2 , 



t=o 2 7 



16817 

120 
304 

'T5" 7 ' 



(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 
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where 7 = K . The first time derivatives are calculated from the boundary condi 
tions 

A 



t=0 



whereas /o is calculated from eq. (5.11) evaluated at t — and fi 



t=o 



t=o 



(2.16) 
is calculated 



from eq. (5.12) using the former result for fo, i.e. 



fo 

A 



t=0 



t=0 



c\H + 2- 



G 



(i)' 



K 



A K\A/2. 



Finally, differentiating eq. (5.11) and evaluating at t = we have 

fo 



2c 2 K 



t=0 



K K V A / 2. 



(2.17) 
(2.18) 



(2.19) 
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